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ABSTRACT 

We use the two-zone model of Cooper & Narayan to study the onset and time evolution of hydrogen-triggered 
type I X-ray bursts on accreting neutron stars. At the lowest accretion rates, thermally unstable hydrogen burn- 
ing ignites helium as well and produces a mixed hydrogen and helium burst. For somewhat higher accretion 
rates, thermally unstable hydrogen burning does not ignite helium and thus triggers only a weak hydrogen 
flash. For our choice of model parameters, these weak hydrogen flashes occur for 10"^ < M/Mfidd < 3 X 10"^ 
The peak luminosities of weak hydrogen flashes are typically much lower than the accretion luminosity. These 
results are in accord with previous theoretical work. We find that a series of weak hydrogen flashes generates a 
massive layer of helium that eventually ignites in an energetic pure helium flash. Although previously conjec- 
tured, this is the first time such bursting behavior has been actually demonstrated in a theoretical model. For yet 
higher accretion rates, hydrogen burning is thermally stable and thus steadily generates a layer of helium that 
ultimately ignites in a pure helium flash. We find that, for a narrow range of accretion rates between the mixed 
hydrogen and helium burst and weak hydrogen flash regimes, unstable hydrogen burning ignites helium only 
after a short series of weak hydrogen flashes has generated a sufficiently deep layer of helium. These bursts 
have fluences that are intermediate between those of normal mixed hydrogen and helium bursts and energetic 
pure helium flashes. 

Subject headings: dense matter — nuclear reactions, nucleosynthesis, abundances — stars; neutron — X-rays: 
binaries — X-rays: bursts 



1. INTRODUCTION 



Type I X-ray bursts are thermonuclear explosions that oc- 
cur on the s urfaces of accreting neutron stars in low-mass X- 
ray binaries dBabushkina et al .il975l; lGrindlav & Heis g[T97l 
Grindlav et al. 1976; Belian et alj 119761: IW ooslev & Taam| 
1976tiJosslil977EiMaraschi & Cavalierelll977l:lLamb & Lamt 



19771 Il978h . They are triggered b y un st able hydrogen or 
helium burning (for reviews, see tewin elai] fT99lTT99l 
ICummingl 120041: IStrohmayer & Bildsten .2001 ). Theoretical 
studies of the onset of type I X-ray bursts suggest that there 
are in general three different bursting regimes separated by 
accretion rate (Fujimoto et al. 1981; Fushiki & Lamb 1987; 
iNarayan & Hey]|l2003l) . At the highest accretion rates, he- 
lium ignites in a hydrogen-rich layer and triggers a mixed 
hydrogen and helium burst. At somewhat lower accretion 
rates, nuclear burning via the hot CNO cycle depletes hy- 
drogen before helium ignites, and so thermally unstable he- 
lium burning triggers a pure helium flash. At the lowest ac- 
cretion rates, hydrog en burni ng itself is thermally u nstable 
(ICzerny & Jaroszynsk i 1980; Ergma & Tufiikovll980l) . 

The vast majority of the type I X-ray bursts that have been 
observed so far are mixed hydrogen and h elium bursts from 
system s accreting at rela ti vely hi gh rates (Ivan Paradijs et alJ 



with the pure helium flash regime. The Wide Field Cameras 
onboard the BeppoSAX satellite observed nine systems that 
exhibited type I X-ray bursts, but from which no persistent 
emission was detected (in 't Zand et al. 1998a; Cocch fet alJ 
119991 1200 It iKaptein et al.ll2000t ICornehsse et alJ i20021. The 
lack of persistent emission suggests that these "burst-only" 
sources have accretion rates M < O.OlMEdd, where Medd 
denotes the mass accretion rate at which the accretion 
luminosity is equal to the Eddington limit (for a review, see 
ICornelisse et al.ll2004l) . The low accretion rates imply that 
the type I X-ray bursts from these sources likely occur in the 
thermally unstable hydrogen burning regime. 

Although it is well understood that hydrogen burning is 
thermally unstable on neutron stars that accrete at low rates, 
it is not immediately apparent what will happen after hy- 
drogen ignites. While thermally unstable hydrogen burning 
will usually ignite helium as well and produce a mixed hy- 
drogen and helium burst, previous investigations have found 
that there exists a range of accretion rates in which ther- 
mally unstable hydrogen burning does not trigger unstable 



helium burni ng and therefore produ c es only a weak hy- 
drogen flash (lErgrna_& Tutukovll 19801: iFujimotoetalJ 1 1981 
iFushiki & Lamblll987l: iPeng et alJ l2007^. The low peak X- 



T988t ICornelisse et al.1 llOOl " iRemiUard et alJ 120061: 
Gallowav et alJ I2006h . This is perhaps to be expected. 



for the predicted type I X-ray burst rate is in general an 
increasing function of the accretion rate, and theoretical 
models predict that the range of accretion rates in which 
mixed hydrogen and helium bursts oc cur is rather large. Re- 
cently, iGalloway & CummingI (l2006l) observed pure helium 
bursts from th e accretion-powered mi llisecond pulsar SAX 
J1808.4-3658 din 't Zand et al.lll998bl) . These are the first 
observed type I X-ray bursts to be unambiguously associated 

Electronic address: |rcooper@cfa.harvard.edu, marayan@cfa.harvard.edu | 



ray luminosity of the weak h ydroge n flashes co uld p ossi- 
bly render them undetectable. iFushiki & LambI (Il987l) and 
iPeng et all d2007i) hypothesize that, within this range of ac- 
cretion rates, a series of weak hydrogen flashes will produce a 
sizable layer of nearly pure helium that will eventually trigger 
an energetic helium flash. 

In this investigation, we use a suitably modifie d version 
of the two-zone burst model of ICooper & NaravanI d2006) to 
study hydrogen-triggered type I X-ray bursts at low accretion 
rates. We find that, for a certain range of accretion rates, ther- 
mally unstable hydrogen burning does not ignite helium, in 
agreement with previous investigations. More importantly. 
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we find that, within the aforementioned range of accretion 
rates, a series of weak hydrogen flashes indeed generates a 
layer of helium that eventually ignites in an energetic helium 
flash, thereby c onfirmi ng the hypothesis of iFushiM & Lama 
(Il987h and Pen g et alJ ^07). We begin with a description 
of the modifications to the two-zone model in ^ We discuss 
the equilibria of the governing set of differential equations and 
analyze their stability as a function of accretion rate in ^ In 
^we integrate the governing equations to determine the na- 
ture of type I X-ray bursts at different accretion rates. We then 
summarize our results and conclude in ^ 

2. THE MODEL 

We assume that matter accretes spherically onto a neutron 
star of gravitational mass M = I AMq and areal radius R = 
10 km at an accretion rate per unit area E as measured in the 
local frame of the accreted plasma. We consider all physical 
quantities to be functions of the column depth S, which we 
define as the rest mass of the accreted matter as measured 
from the stellar surface divided by AttR-. We describe the 
composition of the matter by the hydrogen mass fraction X, 
helium mass fraction Y, and heavy element fraction Z = 1 - 
X-Y. The mass fractions at E = 0, Xq, Yq, and Zq, are those of 
the accreted plasma. In this work, we assume that all heavy 
elements are CNO and that the composition of the accreted 
matter is that of the Sun: Xq = 0.7, Yq = 0.28, and Zq = 0.02. 

To investigate the stability and time evolution of 
the accreted plasma, w e use the two-zone model of 
ICooper & Naravaiil (l2006h which consists of (1) a zone that 
begins at the surface of the star, where E = 0, and extends 
to the depth Eh at which hydrogen is depleted via nuclear 
burning, and (2) a zone that begins at Eh and extends to the 
depth Ehb at which helium is depleted via nuclear burning. 
Hydrogen, helium, and CNO are all present in zone 1, while 
only helium and CNO are present in zone 2. Furthermore, 
we include both hydrogen and helium burning in zone 1, 
but we include only helium burning in zone 2. Tu and 
Zh denote the temperature and CNO mass fraction at Eh, 
respectively, and Thb denotes the temperature at Ehb- For the 
present study we introduce three modifications to the original 
two-zone model. (1) For the low accretion rates considered 
in this work, the temperatures in the accreted layer are often 
below 8 X 10^ K. At these low temperatures, '■'N /3-decay 
competes favorably with '^N(/?,7)'"^0, and so hydrogen burns 
predominantly via the cold CNO cycle (' Wiescher et al.|[T999h 
i2C(/:',7)i^^N(/3+i/)i3C(/:',7)i4N(p,7)i50(/3+z/)i5N(/7,a)i2C, 
the rate of which is determined primarily b y the slow 
tempe rature-dependent reaction '"^N(/?,7)'''0 (.Arnould et al.l 
[T992h . Therefore, we use the more general prescription 
for the hydrogen nuclear energy generation rate en de- 
scribed in Naravan & Hevl (2003), which includes cold 
CNO cycle hydrogen burning, with the exceptions that the 
'^N(p,7) ''*0 and ''*N(p,7)^^0 reaction rates are updated to 
fliose of ICaughlan & Fowled (Il988l) . (2) In this work, the 
outward heat flux due to non-equilibrium electron captures, 
neutron emissions, and pycnonuclear re actions that occur 
in th e both the outer and inner crust (Haens el & Zdunifl 
119901 I2003t iBrownl l2000t iGupta et al.. 2006,) largely sets 
the thermal structure of the acc reted layer (iFujimoto et al.l 
Il984t iHanawa & Fuiimotol [T98?). This is in contrast to the 
accreted layers on neutron stars that accrete at high rates, 
for which mixed hydrogen and heliur n burning set s the 
thermal structure. Therefore, we follow iBrownl (120001) and 
set the outward flux at the base of the accreted layer due to 



deep crustal heating Fbase = 0.1 MeV(E/mu), where ;«„ is the 
atomic mass unit. (3) In their derivation of the two differential 
equations that govern the time evolution of the temperature 
The a nd column depth Ehb o f zone 2, the helium-burning 
zone, ICooper & Naravaiil ^2006) assumed that instantaneous 
changes of the physical quantities Eh, Zh, and Tu of zone 
1 affect the physical quantities Ene and The of zone 2. This 
assumption was physically well-motivated for their work 
because the thermal and accretion timescales of zone 1 
were on the order of the thermal and accretion timescales of 
zone 2, respectively. As we show in ^ however, the two 
timescales of zone 2 in this work are much greater than their 
counterparts in zone 1, so the previously stated assumption 
is no longer valid. Therefore, we drop this assumption when 
we derive the differential equations that govern the current 
model. The result of the new derivation is equivalent to 
simply omitting the dJ^u/dt, dZn/dt, and dTa/dt terms 
that appear in equations (32) and (33) of ICooper & NaravanI 
(2006). The fundamental equations of our two-zone model 
are thus 



dt 
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'^O^'h 



dt 

dTu _ 5 
dt 4C„ 
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4AE 



dY, 



He 



dt 



(5) 

See ICooper & Naravaiil (l2006l) for the definitions of the vari- 
ous symbols. 

We emphasize that, in an effort to construct a model that 
is as simple as possible and yet contains all of the neces- 
sary physics needed to describe hydrogen-triggered type I X- 
ray bursts, we make several assumptions that are not always 
valid. First, we follow ICooper & NaravanI (120061) and assume 
an ideal gas equation of state, but the electrons are degenerate 
and set the equation of state at large column depths. In the 
governing equations ([T]|5]l however, the equation of state of 
the matter shows up only in the specific heat at constant pres- 
sure Cp = 5ks/2mp. The ion specific heat should dominate 
that of the electrons for the range of column depths we con- 
sider, so our approximation is reasonable. Second, we ignore 
the density dependence of the helium nuclear energy gener- 
ation rate ene, and so the rate is somewhat overestimated in 
zone 1 and underestimated in zone 2. However, considering 
the other approximations we have made it should be accu- 
rate enough for our purposes. Third, we set the opacity to the 
constant value k = 0. 1 36 cm^ g"' . This approximation is prob- 
ably the most severe. Electron conduction sets the opacity at 
the largest column depths we consider in this work, and so 
the true opacity can be notably smaller than we assume. We 
discuss the consequences of this approximation in ^ 
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Fig. 1 . — Equilibrium values of the four variables Th, Tro Sh, and Sne as 
a function of /acc = S/Sndd- Zj^'' ~ ^0 the entire range of /acc considered 
in this work. S^'' reaches a local minimum at roughly the same /^cc at which 
is great enough to initiate hot CNO cycle hydrogen burning. 

3. EQUILIBRIA AND THEIR STABILITY 

In this section, we solve for the equihbrium solutions of 
the two-zone model and determine their stability. To find the 
equilibria, we set d/dt = Q,so that equations ([T]|5J become the 



following set of five coupled algebraic equations: 

eH(ZH,7k)I]H = SXo£;^, (6) 

eHe(7k)SH = S(ZH-Zo)£j^e: (7) 

^k-^ke= [eH(ZH,7k) + eHe(7k)]SH, (8) 

eHe(rHe)AI]=S(l-ZH)£;5e, (9) 

^He-^base = eHe(7He)AS. (10) 

We solve equations (l6lfT0]i for the five equilibrium values , 
Z^, r^'', E^'J., and T^'^. When we do this, we obtain 

C = 'f'base+S[^()£H + (l-Zo)£He] , (H) 

Fl5^ = flase+S(l-Z^V*,. (12) 



Thus the flux emitted from the stellar surface, F^, equals the 
flux released via steady-state nuclear burning of all the fuel 
plus the flux fbase due to deep crustal heating, and the flux 
entering zone 1, F^, equals the flux released via steady-state 
nuclear burning of the helium within zone 2 plus that due to 
deep crustal heating. See Figure[T]for plots of T^'^, T^^, , 
and E^l. as a function of the Eddington-scaled accretion rate 

L,, = ^, (13) 

where Sndd = 1.0 x 10^ gcm"^s"'. 

The trends seen in Figure [T] may be understood as fol- 
lows. Hydrogen and helium burning generates energy 
and heats the accreted layer In equilibrium, hydrogen 
and helium burn at a rate oc l^^c, and so the equilibrium 
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Fig. 2. — Upper panel: real parts of the positive eigenvalues as a function 
of /acc. Lower panel: normalized components of the eigenvector correspond- 
ing to the eigenvalue with the greatest real part. The Zh and Shc components 
are negligible and thus are not shown. Both the eigenvalues and eigenvectors 
show that perturbations in 7h trigger bursts for log(/acc) < — 2.5, and pertur- 
bations in Thc trigger bursts for log(/acc) ^ —2.5. 

temperatures T^'^ and Tj^^ are monotonically increasing 
functions of Zacc- For log(/acc) ^ -2.5, the rates at which 
both hydrogen and helium burn are very temperature- 
sensitive, and so the depths at which hydrogen and helium 
deplete, and S^^, are decreasing functions of ^cc- 
For Tj^"* > 8 X 10^ K, however, the proton capture rate 
onto '"'N(ri/2 = 598 s) exceeds that of its /9-decay, and 
so hydrogen burning proceeds via the hot CNO cycle 
'2C(/?,7)i%(p,7)''*0(/3+j/)i'*N(p,7)i50(/3+i/)i5N(p,a)i2c, 
the rate of which is determined by the sl ow /3-decays of 
1^1/2 = 70.6 s) and '^0(fi/2 = 122 s) (iHovle & Fowled 
IT965I) . In this case, the hydrogen burning rate depends 
only on the CNO mass fraction Zh, which is w Zo because 
^He(T^'^) / E^^, the rate at which helium burns to carbon, is 
negligible. Equation ^ thus implies that S^'' c>c S. This 
regime is reflected in the positive slope of S]^(/acc) in Figure 

m 

After we have solved for the equilibrium solution, we con- 
duct a linear stability analysis by assumi ng that all quantities 
vary with time as exp(Af) as described in ICooper & NaravanI 
(2006). The top panel of Figure |2] shows the spectrum of the 
positive real parts of the eigenvalues A, and the bottom panel 
of Figure|2]shows the normalized squared moduli of the com- 
ponents of the eigenvector v corresponding to the eigenvalue 
with the greatest real part, both as a function of Zacc- 

One can ascertain the physical meanings of the eigenvalue 
spectrum by considering perturbations of proper subsets of 
the system of equations ([T]|5]). A linear stability analysis of 
only equations ([T]l and ([3]), which describe perturbations of 
the hydrogen-burning shell, produces two eigenvalues that are 
virtually indistinguishable from the two largest eigenvalues 
of the full system. Similarly, a linear stability analysis of 
only equations (HI and (|5]), which describe perturbations of the 
helium-burning shell, produces two eigenvalues that are virtu- 
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ally indistinguishable from the two lower eigenvalues (which 
vary nearly linearly with /acc in Figure|2]l. Thus the two largest 
eigenvalues that abruptly become negative at log(Zacc) ~ -2.5 
describe the stability of hydrogen burning in zone 1 , and the 
other two nearly linear eigenvalues describe the stability of 
helium burning in zone 2. The clear distinction between the 
two sets of eigenvalues implies that, according to the linear 
stability analysis, there is little interaction between zone 1 and 
zone 2. Physically, this is because E^"* ^ E^^., i.e. the two 
burning fronts are well-separated. One can deduce this also 
from the 5x5 Jacobian matrix of the full linear stability anal- 
ysis. The Jacobian is nearly a block diagonal matrix, with a 
3x3 block matrix in the upper left corner that corresponds to 
the Jacobian of equations ([T]|3]l and a 2 x 2 block matrix in the 
lower right corner that corresponds to the Jacobian of equa- 
tions (|4]|5]l. Although the two off-diagonal block matrices are 
nonzero in general, their matrix elements are negligible rela- 
tive to those of the two block diagonal matrix elements, which 
suggests that the full system is reducible to a set of two essen- 
tially independent systems. 

Only the mode corresponding to the eigenvalue with the 
greatest real part det ermines the behavior of the syste m near 
the equilibrium (e.g.. lGuckenheimer & Holmeslll983h . From 
a quick examination of the largest eigenvalue and its corre- 
sponding eigenvector shown in Figure |2] our linear stability 
analysis suggests that thermally unstable hydrogen burning in 
zone 1 triggers bursts for log(Zacc) ^ -2.5 and thermally un- 
stable helium burning in zone 2 triggers bursts for log(Zacc) ^ 
-2.5. This agree s well with the resul t s of previous theoretical 
investigations (Fuii moto et aLlflOSTl iFushiki & LambI 119871: 
[Narayan & Heyl 2003t). Note that the value of /acc at which 
this transition occurs is in general a decreasing function of 
-Fbase/S, which is assumed to be (0.1 MeV)/OTu in this work. 

Figure [U shows that 7;^'' > 8 x 10^ K for log(Zacc) > -2.5, 
and so hydrogen burns via the temperature-independent hot 
CNO cycle. The hydrogen-burning layer is thus stable to ther- 
mal perturbations, and so the eigenvalues corresponding to the 
hydrogen-burning layer become negative. The thermal eigen- 
value of the hydrogen-burning layer passes through zero at 
approximately the same las^cS-t which reaches a local min- 
imum (e.g.. ,Paczvriski 1983h . 

4. NUMERICAL INTEGRATIONS OF THE TWO-ZONE MODEL 

The linear stability analyses presented in ^describe the 
behavior of the dynamical systems near their respective equi- 
libria. However, the time evolution of a system that is dy- 
namically unstable need not come anywhere close to its equi- 
librium. This is particularly important for systems in which 
hydrogen burning is thermally unstable. In this regime, a 
thermonuclear instability ensues when the column depth of 



the accreted layer E « Ej^ <C 



^He' 



and so the behavior of 



Ehc cannot necessarily be deduced from the stability analy- 
sis of small perturbations around the equilibrium. Moreover, 
the linear stability analysis suggests that, in this regime, he- 
lium burning in zone 2 is unstable as well, which makes the 
true time-dependent behavior of the systems especially diffi- 
cult to predict from the stability analysis alone. In this sec- 
tion, we remedy this by numerically integrating equations ^ 
Is) to study the nonlinear development of the instability and 
the limit cycle behavior of type I X-ray bursts at different ac- 
cretion rates. 

We initiate all numerical integrations by perturbing the 
equilibrium solution, and we continue the integration until the 
system reaches a limit cycle. The limit cycle behavior is com- 
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Fig. 3. — Time evolution of type I X-ray bursts for log(/acc) = —3.0. The 
top panel shows the light curve, and the bottom three panels show the time 
evolution of the five physical quantities. Note that Thc > Tk and Shc > Sh- 
Sudden drops in both Eh and Ehc during bursts indicate that both hydrogen 
and hehum have ignited and thus have triggered a mixed hydrogen and helium 
burst. 

pletely insensitive to the initial conditions. Figure [3] shows 
the limit cycle behavior for a system with log(Zacc) = -3.0. 
At the lowest accretion rates we consider in this work, un- 
stable hydrogen burning via the cold CNO cycle raises the 
temperature of the accreted layer to such an extent as to trig- 
ger unstable helium burning as well. The sudden drops in 
both Eh and Ehc that correspond with sudden increases in 
the temperatures Tu and and the outward flux Fu in- 
dicate that both hydrogen and helium ignite in every burst. 
These bursts are the well-known mixed hydrogen and helium 
type I X-ray bursts triggered by ther r nally unstable hydro- 
gen burning dErgma & Tutukovlll980t Fuiimoto et al J [19811 : 



IFushiki & Lambiri987l:lPeng et alJl2007l)" 



For somewhat higher accretion rates, theory predicts that 
the increase in temperature due to thermally unstable hydro- 
gen burning is insufficient to ignite helium. Previous authors 
have speculated that a series of weak hydrogen flashes allows 
for the accumulation of a large layer of helium and eventually 
leads to a strong helium flash (lTaamlll9"8Tt IFushiki & LambI 
1987; Peng et al. 2007). We present in Figure S] the results 
of a calculation for a system with log(/acc) = -2.7 that con- 
firms this suggestion: if thermally unstable hydrogen burning 
does not ignite helium, helium will accumulate up to a column 
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Fig. 4. — Time evolution of bursts for log(/acc) = -2.7. Tlie top panel on tlie left shows the light curve, and the bottom three panels show the time evolution of 
the five physical quantities. Note that 7hc > Th and Ehc > Sh- The increase in Th due to thermally unstable hydrogen burning is insufficient to trigger helium 
ignition. Therefore, Shc increases during a long series of weak hydrogen flashes until helium ignites and produces an energetic pure helium flash. The set of 
panels on the right illustrates the time evolution of the system just before and after an energetic pure helium flash. 



depth Shc ^ 10^^ g cm"^ and eventually ignite in an energetic 
helium flash. Subsequent weak hydrogen flashes ensue until 
a sufficient amount of helium can accumulate and once again 
trigger an energetic helium flash. To our knowledge, this is 
the first study in which this limit cycle behavior has been ex- 
plicitly demonstrated. 

The bursting behavior of systems with accretion rates that 
lie in a small range between those of the two regimes de- 
scribed above can be somewhat chaotic. Figure |5] shows the 
light curve Fh and the time evolution of the helium column 
depth She for five different accretion rates that span the transi- 
tion region between the two aforementioned bursting regimes. 
As we increase the accretion rate from log(Zacc) = -2.93 to 
-2.87, a comparison between and Eee illustrates that the 
bursting behavior evolves from that of exclusively mixed hy- 
drogen and helium bursts, to a combination of both mixed 
bursts and weak hydrogen flashes, to predominantly weak hy- 
drogen flashes with infrequent but unusually energetic mixed 
bursts. The left panel of Figure |5] shows that the amount of 
helium fuel burned in the mixed bursts increases with accre- 
tion rate. Thus the burst fluences of the mixed bursts in this 
transition region are intermediate between the mixed bursts 
at slightly lower l^^c and the energetic pure helium bursts and 



slightly higher Zacc- 

For log(Zacc) ^ -2.5, Figure |2] implies that hydrogen burn- 
ing is thermally stable. Consequently, steady-state hydro- 
gen burning generates a sizable column of nearly pure he- 
lium as accretion ensues until the helium at the base of 
the accreted layer ignites. Figure |6] depicts the limit cy- 
cle behavior of such a system with log(Zacc) = -2.0. The 
existence and nature of these pure he lium-shell flashes are 
well-established in theoretical models ('Fuiimoto et al." 19811 
Hanawa & Fuiimoto 1982; Fushiki & Lamb 1987; Bildste^ 
1998; Gumming & Bildsten 2000; Narava n & Heyl 120031 
Woosley et al. 2004; Gumming et al. 2006). Note fliat the col- 
umn depths Shc at which helium ignites sho wn in Figures ffl 
and|6]are g enerally lower than those found by iGumming et alJ 
(2006) and iPeng et al.l (l2007h . This is due to our simplifica- 
tion that the opacity k = 0.136 cm^g"' throughout the entire 
layer. At these large column depths, electron conduction sets 
the opacity, and so the true k is lower than our assumed value. 
Since the opacity is lower, a larger column of accreted matter 
must be accreted before helium burning can become thermally 
unstable. 

We plot the recurrence time, burst energy, a-value, and 
maximum Eddington-scaled flux of type I X-ray bursts as a 
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Fig. 5.— Burst light curves (left) and Shc (right) for log(/acc) = -2.93, -2.92, -2.89, -2.88, and -2.87, from top to bottom. For log(/acc) = -2.93, the Th 
increase due to thermally unstable hydrogen burning is always sufficient to ignite all of the helium. For slightly larger values of l^c , unstable hydrogen burning 
triggers complete helium ignition only if Shc is sufficiently large, and the ensuing mixed hydrogen and helium bursts are therefore less frequent and more 
energetic. For log(/acc) = -2.87, unstable hydrogen burning never ignites helium. 



function of l^^c in Figure |7] All quantities are as measured by 
an observer at infinity. For -2.9 < log(Zacc) ^ -2.5, both en- 
ergetic pure helium flashes and weak hydrogen flashes occur 
In this regime, the solid lines in the four panels denote the 
burst properties of the helium flashes, and the dot-and-dashed 
lines denote the burst properties of the hydrogen flashes. The 
recurrence time is defined as the time between two successive 
flashes. We calculate the burst energy by determining the total 
amount of hydrogen and helium fuel that exists just prior to 
a burst and assuming that all of the fuel burns during a burst. 
Since essentially no helium burns during a weak hydrogen 
flash, however, we determine the burst energy of such a flash 
by instead assuming that the burst consumes only hydrogen. 
The parameter a is defined as the accretion energy released 
between successive bursts divided by the nuclear energy re- 
leased during a burst. For the range of accretion rates in which 
both energetic pure helium flashes and weak hydrogen flashes 
occur, the larger a in Figure |2]is the value an observer would 
measure if only the energetic helium flashes were detected, 
and the smaller a is the value an observer would measure if 
the weak hydrogen flashes were detected as well. Since the 



peak flux of weak hydrogen flashes is usually much less than 
the accretion luminosity, as is evident in the bottom panel of 
Figure |2l it is possible that weak hydrogen flashes would be 
undetected. Note that, as Figured illustrates, the properties 
of the weak hydrogen flashes for a given l^^c vary somewhat 
with time. Consequently, the quantities plotted in Figure [T] 
that describe weak hydrogen bursts are only approximate. 

The ranges of accretion rates within which the various 
bursting regimes lie depend upon the value of the outward 
flux at the base of the accreted layer Fbase, which we set to 
0. 1 Me V(S / mu). However, for the low accretion rates we con- 
sider in this work, Fbase may be significantly larger (Brow^ 
12004 . When we instead set Fbase = 1.0MeV(E/mu), we find 
that weak hydrogen flashes occur for -3.95 < log(/acc) ^ 
-2.55. This range is lower that that shown in Figure|7] and it is 
qualitatively consistent with the results of lPeng et al.l (l2007h . 

Our numerical integrations expose a fundamental weakness 
of linear stability analyses. While one can accurately deter- 
mine whether or not a system is stable to perturbations using 
a linear stability analysis, one cannot necessarily infer with 
any confidence the nature of a system that is unstable. For 
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Fig. 6. — Time evolution of bursts for log(/acc) = —2.0. The top panel shows 
the light curve, and the bottom three panels show the time evolution of the 
five physical quantities. 7h ^ 10** K, so hydrogen burns via the temperature- 
independent and thus thermally stable hot CNO cycle. Consequently, stable 
hydrogen burning generates a large and nearly pure helium layer that eventu- 
ally ignites in an energetic helium flash. 

example, numerical integrations demonstrate that a transition 
between hydrogen-triggered mixed bursts and weak hydrogen 
flashes occurs at log(Zacc) ~ -2.9, but neither the eigenvalues 
nor the eigenvector shown in Figure |2] indicate any change in 
the physical behavior of the system near this accretion rate. 
This weakness exist s even in the mor e complete global linear 
stability analysis of .Naravan & HevI (12003) . In that work, the 
authors found that hydrogen burning is thermally unstable at 
low accretion rates as well, but again they were unable to de- 
duce the nature of the instability from their results. Thus, one 
must use caution when inferring the time-dependent behav- 
ior of a system that is unstable according to a linear stabiUty 
analysis. 

5. SUMMARY AND CONCLUSIONS 

Using a suitably modified version of the two-zone model 
of ICooper & NarayanI (l2006h . we have studied both the on- 
set and time evolution of hydrogen-triggered type I X-ray 
bursts on accreting neutron stars. At the lowest accretion 
rates we have considered, thermally unstable hydrogen burn- 
ing triggers thermally unstable helium burning and thus pro- 
duces a mixed hydrogen and helium burst. For somewhat 
higher accretion rates, thermally unstable hydrogen burning 
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Fig. 7. — Recurrence time, burst energy, ratio of accretion to burst energy 
a, and peak flux of type I X-ray bursts as a function of /acc, from top to 
bottom. For —2.9 < log(/acc) ^ -2.5, both energetic pure helium flashes and 
weak hydrogen flashes occur. In this case, the sohd line denotes the burst 
properties of the helium flashes and the dot-and-dashed line denotes the burst 
properties of the hydrogen flashes. The chaotic behavior discussed in 0and 
Fig. [5]causes the sudden jumps in the panels that occur in the range -2.93 < 
log('acc) < -2.87. The dashed line in the bottom panel denotes the luminosity 
due to accretion. Note that the peak flux during most weak hydrogen flashes 
is much less than the accretion luminosity. 



does not ignite helium and thus triggers only a weak hydro- 
gen flash, in agreement with previous studies. For our choice 
of model parameters, these weak hydrogen flashes occur for 
10"^ < M/M fidd 5. 3 X 10"l We find fliat, in accord with th e 
predictions of iFushiki & LambI (Il987h and lPeng et all (l2007h . 
subsequent weak hydrogen flashes generate a sizable layer 
of nearly pure helium that eventually ignites in an energetic 
pure helium flash. This is the first time this bursting behavior 
has been seen in theoretical models. For yet higher accretion 
rates, hydrogen burning is thermally stable and thus steadily 
generates a layer of helium that ultimately ignites in a pure he- 
lium flash. In addition, we find that there exists a small range 
of accretion rates near the boundary between the mixed hydro- 
gen and helium burst and weak hydrogen flash regimes where 
thermally unstable hydrogen burning ignites helium only after 
Shc is sufficiently large. The resulting mixed hydrogen and 
helium type I X-ray bursts have fluences that are intermediate 
between the normal mixed hydrogen and hehum bursts and 
the energetic pure helium flashes. 
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The simplicity of our two-zone model has enabled us to 
survey a wide range of accretion rates and to illustrate the ba- 
sic physics of hydrogen-triggered type I X-ray bursts in the 
various bursting regimes. However, there are several pieces 
of physics we have not considered that certainly would affect 
the nature of type I X-ray bursts at the low accretion rates 
considered in this work. First, the timescale for sedimenta- 
tion and diffusion of heavy elements in the accreted layer is 
often less than the accretion timescale (Wallace et al. 1982; 
iBildsten et all 119921: iPen^ et al.l l2007h . Peng et al. (2007i) 
demonstrate that the sedimentation of CNO ions can have a 
notable effect on the ignition conditions of type I X-ray bursts 
at very low accretion rates. Second, our effective nuclear 
reaction network consists only of hydrogen burning via the 
cold and hot CNO cycles and helium burning via the triple-a 
process. During a mixed hydrogen and helium burst, how- 
ever, hydrogen will burn p redominantly via the rp-process of 
[Wallace & WooslevI (11981 ). Consequently, the light curves of 
these bursts would be somewhat different than those shown in 
Figure|3]if rp-process hydrogen burning were included. 

A potentially serious issue that we are currently unable to 
address, but that could have significant effects on the ignition 
of type I X-ray bursts at low accretion rates, is the destruction 
of CN O ions via nuclear spallation reactions ( IBildsten et al.l 
|1992|) . The accretion disks around neutron stars that accrete 
at low rates are Ukely truncated (e.g. ,Done.2002;,Barret 2 004). 



Within the truncated disk, the accretion flow is op tically thin 
and quasi-spherical (see, e.g. iNaravan & Yilll995h . and thus 
the accreting ions have large radial velocities when they reach 
the neutron star surface. Observationally, this radiatively 
inefficient accretion flow occurs in neutron star low-mass 
X-ray binaries when M < O.OlMEdd (Mendez etaH 119971: 
iGierlinski & Donel2002h . the range of accretion rates appUca- 
ble to this work. The directed kinetic energy of the accreting 
hydrogen and helium may be sufficient to destroy nearly all 
of the ions in the envelope that are heavier than helium. The 
degree to which spallation would affect type I X-ray burst ig- 
nition is unclear Compositional inertia of the neutron star en- 
velope could potentially counteract any d epletion of CNO due 
to spallation (Avasli & Joss 1982; Woos lev & Weaved [19841; 
iTaam et alTi993ciWooslev et al.ll200 4). However, if the heavy 
elements produced in a given burst sediment out of the enve- 
lope before the next burst is triggered, compositional inertia 
may be negligible, in which case spallation again becomes 
important. This issue in particular is worthy of further inves- 
tigation. 
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